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The light-cone Hamiltonian, incorporating the nonpertur- 
bative dynamics of the qq system connected by the string is 
solved numerically. The spectrum is shown to coincide with 
that of the center-of-mass Hamiltonian within the accuracy 
of computation, displaying the expected degeneracies of the 
string states, but an overall shift due to different treatment 
of Z-graphs is obtained. The nonperturbative wave functions 
are calculated directly from the light-cone Hamiltonian for 
the first time, thus allowing the explicit estimation of non- 
perturbative effects in the parton model language. In this 
way one obtains the string contribution to the parton distri- 
butions and in the formfactors and structure functions. 



I. INTRODUCTION 

The light-cone description of the hadron wave func- 
tions is widely used, since it allows to get a direct con- 
nection to the parton model and its QCD improvements 
0. The latter however are mostly perturbative and 
nonperturbative contributions are introduced via OPE 
and QCD sum-rules. In this way the concept of the 
string - the main nonperturbative QCD phenomenon - 
is totally lost. On physical grounds it seems that the 
string is the essential ingredient of the dynamics for large 
(r > 0.5 /to) distances, and it would be interesting to 
understand its contribution to the light-cone wave func- 
tions, formfactors, structure functions etc. 

In particular, what is the QCD string in the parton 
language? Should one associate it with the gluon con- 
tribution as an assembly of gluons compressed inside the 
string - or with the constituent quark mass? 

There two contrasting points of view have been pro- 
posed already decades ago ||||. In 0j the sea quarks 
and gluons enter as separate entities and one could asso- 
ciate the gluon distribution with the string in the same 
way as photons with the Coulomb field of the charge in 
the Williams- Weizsacker method. In contrast to that in 
||] the quarks have been considered as constituents with 
structure, and string does not appear separately. Re- 
cently a quantitative analysis was performed Q of quark 
distribution in the pion starting from that in the nucleon 
and assuming the same internal structure of quarks in 
nucleon and in pion. It is still an open question how 
the structure of the constituent quarks is formed, and to 
which extent it can be explained by the adjacent piece 



of the string. This problem can be elucidated partly in 
the present approach, since our light- cone Hamiltonian 
contains the string on the light cone explicitly. It allows 
to separate the contribution of the string to the parton 
distribution, in particular to the momentum sum-rule. 
There is another source of structure in the constituent 
quark - chiral symmetry breaking which creates the chi- 
ral mass of quark. This problem will not be discussed 
below, see e.g. ||. 

There are many other questions which can be asked 
from a decent light-cone Hamiltonian, incorporating qq 
system with the proper light-cone nonperturbative dy- 
namics. An important advantage of the light-cone wave 
functions is that they allow to calculate formfactors and 
structure functions directly, without additional boosts, 
which typically are the dynamical ones, i.e. require the 
use of the Hamiltonian. 

Another point to clarify is the comparison of the proper 
light-cone dynamics incorporated in the wave-function 
of the light-cone Hamiltonian with the popular an- 
zatz where one is simply using the cm. wave function 
expressing it through the light-cone variables (jj . 

But the first question to answer is the comparison of 
the cm. and light-cone spectra. In the cm. system the 
Hamiltonian of the spinless qq system, connected by the 
minimal (i.e. without vibrations) string was obtained 
in The backward-in-time motion or so called Z- 
graphs are suppressed there because Z-graphs contains 
backward moving string and the corresponding increase 
in the action damps the amplitude. On the light-cone Z- 
graphs are treated differently and do not appear in lowest 
orders of perturbation theory and comparing spectrum 
with that of cm. one can visualize the difference due to 
the different contribution of Z-graphs. 

Another typical feature of the light-cone Hamilto- 
nian should be mentioned. As we shall see in Sec- 
tion 2 the form found in explicitly contains the z- 
axis component of the angular momentum, L z which is 
not Lorenz-invariant and hence the light-cone spectrum 
should demonstrate degeneracy for different values of L z 
with the same L 2 . This degeneracy will be explicitly 
demonstrated below. Therefore all states of the system 
including angular excitations may be calculated using the 
Hamiltonian with L z = 0. It is known (see J^] and also 
Q) that the problem under consideration possesses the 
dynamical symmetry which implies that quasiclassical 
frequencies of the radial and orbital motions are approx- 
imately proportional to each other with the coefficient 
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2. Thus the light-cone Hamiltonian should also manifest 
this dynamical symmetry. 

The plan of the paper is as follows. In Section 2 we 
start with the light-cone Hamiltonian in 3+1-dimcnsions 
from H and prepare it for numerical computations. As in 
we simplify the matter considerably, neglecting spins, 
perturbative gluon exchanges and additional quark pairs. 
The first two points are unessential at large distances, 
which are of the main interest for us here but seriously 
contribute to the lowest mass mesons. We plan to in- 
clude these effects in next publications. Sea quarks are 
expected to contribute some 10% effect and disappear in 
the large N c limit. As it is, we are confined to the valence 
quark sector of the Fock column wave function. 

In Section 3 general properties of solutions and the 
physical meaning of parameters fj,i , x and p are discussed. 
We also demonstrate the limiting procedure reducing 
the 3+1 Hamiltonian to the 1+1 problem, following the 
method of || and compare properties of 3+1 and 1+1 so- 
lutions. In Section 4 the formfactors and structure func- 
tions are defined through the solutions of the light-cone 
Hamiltonian and physical limiting cases are discussed. In 
Section 5 the numerical procedure used for the solution 
of the eigenvalue equation is explained. In Section 6 we 
present numerical results and discuss them from different 
points of view, e.g. comparing center-of-mass spectrum 
and wave functions with those of the light-cone. Special 
attention is paid to the procedure used in literature 
where the cm. wave function is kinematical expressed in 
terms of the light-cone variables. A short conclusion and 
perspectives are given in Section 7. Two appendices, A 
and B serve to illustrate the derivation of the Hamilto- 
nian (|l|) and the final form of equation ( |l3"| ) to be solved 
numerically. 



II. THE GREEN'S FUNCTION AND THE 
HAMILTONIAN 

We consider the relativistic quark-antiquark pair with 
the masses mi and TO2 connected by the straight-line 
Nambu-Goto string with the string tension a in 3+1 
dimensional space-time, see [0 and Appendix A of the 
present paper for details. We start from the expression 
for the classical light-cone Hamiltonian function in 3+1 
dimensions 
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(1) 



In Appendix A we briefly derive this expression, using 
H and omitting unessential details. Here pi,^2 are ein- 
bein fields, playing the role of P+ momenta of the par- 
ticles, p, = ^i^ 2 /(mi + M2) > v = Jo d(3v((3) and 
is the einbein field with the physical meaning of the 



P+ -momentum, carried by the string. As it has been 
mentioned above, the Hamiltonian explicitly depends on 
L 2 = p\r\ — (p±r±) 2 , the corresponding mass param- 
eter in the denominator in (^) is equal to 



a = Ml (l - x) 2 + fi 2 x 2 + J df3v{!3){f3 - x) 2 (2) 



The variable x in the above expression is defined as (see 
Appendix A for details): 



/xi+ < (3 > v 
P+ 



(3) 



where 

P+ = Mi + M2 + v a (4) 
is the light-cone total momentum of the system. Also 



7 = i/ < £ > 



Mi + M2 



(5) 



where < >= d{3fiv{{i)l df3v(p). 

The corresponding Green's function is defined as an in- 
tegral over the dynamical fields as well as over the cinbein 
fields: 



G(xx;yy) — J D^iD^DvDR^Dr^e 



(G) 



where A is the euclidean action. Separating out the cen- 
ter of mass motion one defines the Hamiltonian from the 
corresponding Minkowskian action A M : 



A 



dz + L M , H = p ± r ± - L 



A I 



(7) 



withpj^ = dL M /df±. This is the light-cone Hamiltonian 
of our problem ([!]). 

As the next step one should quantize the classical 
Hamiltonian function. Before doing it, one should choose 
the appropriate set of dynamical variables. Three einbein 
fields Mi)M2> v introduced above play different dynamical 
roles. In the nonrelativistic case m 1 ,m 2 >• \fa (and 
therefore for the free particles) the dependence of Hamil- 
tonian (and wave functions) on v can be correctly found 
by minimization procedure with v taken as a classical 
variable. This is in its turn a consequence of the fact, 
that string in our approach is the minimal one - it has 
no internal degrees of freedom and may only stretch or 
rotate as a whole. 

On another hand, fj,i and [12 on the light cone play the 
role of legitimate quantum dynamical degrees of freedom 
and can be expressed through x and P+ as in (|J) and (||). 

There are two canonically conjugated pairs {p*x,r±y 
and {x, (P + r_)} (see Appendix A). We introduce also 
a new dimensionless variable y instead of v : y(f3) = 
v{f3)/P + . It satisfies the obvious condition: < y < 1. 
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This variable depends on r± as well as on (P + r_). Rigor- 
ously speaking one should extract this dependence by the 
minimization of the Hamiltonian with respect to y as it 
has been explained above. Instead an easier (however ap- 
proximate) way is chosen. On physical grounds it can be 
shown that for small r\ one has yo = J Q yd/3 ~ const ■ r\ 
so that linear string energy density yo/r^_ stays finite if 
r\ — > 0, and using this property one can reproduce the 
correct 1+1 limit, namely 't Hooft equation [{Tcfl , from 
the 3+1 light-cone Hamiltonian (0) (see M and the next 
section). On the other hand, if distances are increasing, 
yo tends to some limiting value which is determined by 
the virial theorem arguments. So one can parameterize y 
introducing several parameters and replace the minimiza- 
tion in the functional sense by the ordinary minimization 
with respect to these parameters. We have chosen the 
simplest 2-parametric form for y: 



V 



yt 



l + at 



(8) 



where t — r\ and y and a are free parameters. The re- 
quirement < y < 1 leads to the restriction < y < a. 
Let us stress again that this parameterization is the mat- 
ter of convenience and all physical results are determined 
from the requirement that every energy level should have 
its own minimum. 

It is easy to see from (^), that the x-variable is the 
part of the total momentum, carried by the first quark 
itself and a part of the string, " belonging" to this quark. 
Rewriting (^) in the form: 

^-=X~y<[3> ; ^ = (l- X )-y(l- < > ) ( 9 ) 
*+ "+ 

one can conclude, that for the given y the variable x may 
vary in the following limits: 



y < (3 > < x < l-y (1- < >) 



(10) 



It is more convenient to make rescaling to a new variable 
p which varies from zero to unity: 



x = y < /3 > +(1 - y)p ; < p < 1 



(11) 



The quantity we are really interested in is the mass op- 
erator squared: 



M — 2 H P+ 



(12) 



and for M 2 we have the Schrodinger equation of the fol- 
lowing form: 

Ai tp" + A 2 ^ + A 3 ip' + A±ip' + A 5 ip + A 6 ip = M 2 ip 

(13) 

here tp = ip(p,r±) and two independent variables are 
p G [0, 1] and t = r\ e [0, oo) ; tp' = dtp/dt; ip = dtp/dp 
; the coefficients Ai depend on p and t. The derivation 



of ( p^| ) from the Hamiltonian (Q) as well as the explicit 
form of the functions A t may be found in Appendix B. 
The following boundary conditions are to be imposed: 

V>(0,t) =V(M) = , V(p,i^oo) = (14) 

with the normalization of the solution: 



dp / d 2 r ± \ 1 p(p 1 f 1 _)\ 2 = 1 



(15) 



We discuss properties of solutions of (|T^) in the next 
section. 



III. GENERAL PROPERTIES OF THE 
SOLUTIONS. 

In this section we shall discuss the physical meaning 
of our variables x, p, p and y and some properties of the 
solution tp(p, t). Rewriting (0) as 



, = Mi M2 



y 



(16) 



one can conclude that Pi/P+ and y have the meaning of 
the parts of the total momenta P+ , associated with the 
i— th quark and string respectively. This identification is 
also supported by the form of the Hamiltonian (|l]), where 
masses enter as m 2 /2/ij so that in the infinite-momentum 
frame it becomes m 2 /2p 1 + . The similar situation occurs 
in the cm. system where p,i and v have the meaning 
of the energy, associated with the quarks and the string. 
One notes the essential difference in the way these ein- 
bein fields are treated in the cm. Hamiltonian and on 
the light-cone. In the former case pi enter on the same 
grounds as v (or y) and are to be found by minimization 
of the Hamiltonian. In contrast to that on the light-cone 
Pi,P2 are expressed through y and x as in (|3j) and x is 
the canonical momentum, conjugated to the coordinate 
P+r_ ( see Appendix A), namely: 



fx, P+r\ = -i 



(17) 



Hence pi are promoted to be quantum operators with 
average values to be found from the solution of the 
"Schrodinger equation" ([l3]) Without the string, when 
y = one has p\ — xP+ ; p,2 = (1 — X )P+ an d the 
correspondence with the naive parton model is trivial. 
However when y =/= there are no a priori reasons to 
associate our x with the Feynman x parameter, but in- 
stead one could call Pi/P+ the Feynman x parameter and 
associate y with the Feynman parameter for the string. 

However the strong interaction between the string and 
quarks taken into account in the light-cone Hamiltonian 
(|l|) invalidates the literal interpretation of our results in 
the free parton language - one can only interpret the 
averaged values of pi,p2 together with the minimized 
value for y. 
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Then one can consider the momentum sum-rules for 
the quark structure functions (%f) where the quark 
Feynman parameter xp normalized to the interval [ 0,1 1 
should be associated with the parameter p from (|ll| ) 
which is connected with pi as follows: 

^- = {i-v) P ; ^- = {i-y){i-p) (18) 

Therefore as we shall discuss later the average momen- 
tum carried by the quarks is equal to (1 — y) while the 
rest, y is carried by gluons assembled into the string - 
nonperturbative gluons. This interpretation of the miss- 
ing momentum (around one half in experiment pj] ) as 
due to the string was suggested already in Now we 
are in the position to calculate this quantity and compare 
to the experiment which we shall do in Section 6. As we 
shall see, the main missing ingredient in our Hamiltonian 
and wave functions is the sea quarks and "sea gluons" - 
i.e. gluons which excite the string and make it vibrating 
instead of keeping the straight-line form we have con- 
sidered. From this point of view our wave function, the 
solution of ( |l3| ) is the first upper entry in the Fock col- 
umn, containing states with 1,2,... gluons, "sitting on the 
string" and exciting it to the first, second, etc. vibrating 
string state. It is the whole Fock column which is re- 
sponsible for the Regge behaviour of cross-sections and 
Regge asymptotic of structure functions near p — > 0. In 
this paper we concentrate on the properties of the sector 
with two valence quarks and minimal (i.e. nonvibrat- 
ing) string, leaving the general discussion of the QCD 
reggeons and the pomeron to a future publication. 

Let us now turn to the properties of the wave functions 
and, in particular, their behaviour at the boundaries. In- 
serting in (13) ip(p) ~ p a and (1 — p) a near the points 
p = and 1 respectively one obtains the only solution 
q = 1 at all t hence the boundary conditions (|l4|) are 
satisfied. Similarly, representing ip(p,t = r\) at large t 
in the form 



i>(j>,t) ~ A(p)exp(-~/(p)t$) 



(19) 



yields £ = 1 and some complicated equation for j(p). 

Finally, we discuss in this section the limit of the (1+1) 
't Hooft equation |l(J . Following the arguments given in 
H we look for the limit t — > with y — yt/{\ + at) — > 
and y fixed and to be found from the minimization. As 
a result, one obtains: 



Hi 



Mi 



P2 



1/09) 



dp + / y(P)dp 



P+r 2 _ 



Mi 



M2 . 

(20) 



Note that the string part in the sum (16) disappears 



because of y = but y is nonzero and ensures the string- 
like potential dynamics. Taking the extremum of ( p3 ) 
with respect to y(P) and introducing x instead of p\, pi 
where p\ = P+x,p2 = P+0- ~ x ) one obtains: 



M 



2P+H 1 



1-x 



2a\P+r-\ 



(21) 



From (pl|) one easily deduces the 't Hooft equation p0| 



dyip(y) 



M 2 tp(x) (22) 



x 1 — x ) T y ' 2n J (x — y) 2 
Thus the 't Hooft solutions pi which approximately are: 



VVi ( x ) = Cn sin wnx 



(23) 



can be thought of as limiting case t — > of our solutions 
but one should have in mind that the extremum value of 
y leading to the 't Hooft equation (|2|) in general does 
not coincide with the extremum value in our equation 
( |l3| ) obtained for all nonzero t. Still the behaviour at 
the boundary points x = 0, 1 (p = 0, 1 in our case) and 
the general pattern of our solution ijj(p, t = 0) is closely 
resembling those of 't Hooft solutions. 



IV. FORMFACTORS AND STRUCTURE 
FUNCTIONS 

The formfactor of the bound system of two spinless 
quarks can be expressed through the wave-function in 
the momentum representation, depending on x and on 
the relative momentum p± as follows (for equal mass 
quarks) O: 



F{<f) = J ip(x,k x )ip*{x,kx + (1 -x)q x )dxd z k x (24) 

Our wave function is defined in the mixed space p,r\= t 
hence one obtains: 

oo 1 

2\ -It !././•_ j.\|2 



Jq <^ qVt 



F( Q r , )=7r j j m P ,t)\ 





1. 

i- p + y{p- p 



| dpdt 



(25) 



where we have exploited relation ( |ll| ) with < P >= 1/2. 
For q 2 = the usual normalization follows F(0) = 1 since 
our wave function is normalized as 



\ip(p,t)\ 2 dpdt = 1 



(26) 



o o 



Note that the effect of the electromagnetic current, 
which interacts only with the charged quark (and anti- 
quark) and does not interact directly with the string, en- 
ters @ through the factor (l-x)q = [l-p+y(p-l/2)]q, 
and not through p or p. This is the result of the fact, that 
the parameter x governs the distribution of the momen- 
tum between the center of mass and relative coordinates 
R± and r ± . 
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Let us notice, that because we neglect perturbative 
gluon exchanges, our formfactor does not have the famil- 
iar quark asymptotic F(q 2 ) ~ const/ q 2 for q 2 — > oo, but 
rather corresponds to the nonperturbative part which is 
believed to dominate the low-g region (actually the ex- 
perimentally accessible region of today) [[L2| . 

From (|2^) and the behaviour of ip(p,t) at large t and 
p — > 1, equation (|l9|), one easily obtains that 



F(q 2 ) 



const 



oo 



(27) 



We turn now to the structure function (or quark dis- 
tribution function). As well as in case of the formfactor, 
since as it is produced by the external current it is natural 
to define it a function of x, but to compare it with stan- 
dard results normalized to the interval [0,1], we choose 
to define it as a function of p. 

Thus, the parton distribution inside our meson wave 
function can be written as: 



q(p)=7T / \i>(p,t)\ 2 dt 



(28) 



and the normalization of q(p) due to (Eq) is 



q(p)dp = 1 . 



(29) 



The momentum sum-rule, using (|16| ) and (|18|) is (for 
equal masses of quark and antiquark) 



q(p)p{l- < V >)dp + / q(p) <y>dp = l (30) 



where < y > is obtained by integration with \ip{p, t)\* 



<y>- 



y\^(p, t)\ 2 dpdt 



(31) 



Comparing with the standard energy-momentum 
sum-rule [jlj containing the gluon distribution g(x), one 
may associate it with y as follows J Q g(x)xdx y > . 

From (||(| [l0|) one can see that both valence quarks 
together carry the part of momentum P + equal to 1— < 
y >, while the string carries on average < y > P+ ■ 

The essential difference of the minimal string from the 
gluon distribution g(x) is that the former is fixed in the 
quasiclassical einbein formalism and has no dispersion - 
unlike free gluons which can carry any part of the total 
momentum P + . 

Let us discuss finally the correspondence between q{p) 
and F(q 2 ) (the Drell-Yan-West duality g§). If F(q 2 ) 
behaves as (q 2 )~ n at large q 2 , then F 2 {x) should decrease 
as (1 — x) p for x — ► 1, where In = p + 1. In our case from 
( p7j ) n — 3/2, hence p — 2 and this agrees with (28), 
since tp(p) ~ (1 — p), and q(p) ~ (1 — p) 2 . 



V. NUMERICAL SOLUTION 

The equation ( |l3| ) is the linear differential equation 
of the elliptic type. Boundary conditions are given in 
(|l4|) and are compatible with the singularities of coeffi- 
cients A4. We have not succeeded in finding any analytic 
function, describing the main features of the solution, 
for example, one can demonstrate that simple ansatz 
ip(p,t) = 4>(p) ■ £(t) is not a solution of (|l3|). Therefore 
we admitted two numerical procedures: a direct solution 
of (|l3| ) using discretization and a variational ansatz with 
the expansion of the solution into a complete set of func- 
tions of two variables. We have found that the former 
procedure turned out to be unstable and not accurate. 
We discuss below only the numerical results obtained in 
the framework of the second method. 

The ^ function in ( fL3| ) depends on parameters mi, TO2 
(bare quark masses), a (string tension), 2 parameters y 
and a defining the string parameter y (see @). The \1/ 
also depends on quantum numbers: L z , total momentum 
L and radial quantum number N r . The eigenvalue of 
( |l3| ) is the bound state mass squared M 2 (L Z , L, N r , y, a). 
We expand ^-function into the complete set of functions 
(k = L z ,L,N r ): 



OO OO 



V K (j), t) = exp{lL z 4,) J2 C ranl^ 



mnl 



(32) 



where 

ipmni(p, t) = sin(7rmp) i 5 L l n (et) exp(-et/2). (33) 

Here <f> is the azimuthal angle of the corresponding an- 
gular motion in the xy-plane and L l n - Laguerre poly- 
nomials. The variational parameter e and the definition 
I = L z were also introduced. The basis functions are 
orthonormalized as: 



OO 2-7T 1 

dt 



000 



dp^>n>l'(P, t )4>rnnl(p,t) 



AAA 77 ( n + 1 )- ..>,, 
Om'm On'n Of l [+1 — (64:) 



2e l+1 nl 
and the total function is normalized as: 

J dV\t>\ 2 = 1 dV = -dpdtd4> 



(35) 



Note, that each function ip mn i satisfies boundary condi- 
tions (|l4|). We have used typically seven sine functions 
(to = 1,..,7) and seven Laguerre polynomials, which 
was compatible with the accessible computer time. The 
eigenvalue equation (13) was transformed into 49 x 49 
matrix equation and solved in a standard way. In some 
cases in order to check sensitivity of the procedure also 
other sets have been used, i.e. 5x9 instead of 7 x 7 and 
consistency of results was confirmed. 
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One important remark should be done here. The ba- 
sis we have used is not large enough to analyze highly 
excited states with quantum numbers much more than 
unity because of loss of the accuracy. But in this region 
powerful quasiclassical methods may be applied in order 
to get the solution. Indeed, one of the main questions we 
are going to solve is to compare the spectra of light-cone 
Hamiltonian and center of mass one for the same dynam- 
ical problem: two massive spinless quarks connected by 
the straight-line string. In 1+1 case it is known for a 
long time |Q , that light-cone and center of mass Hamil- 
tonians for this problem are quasiclassically equivalent. 
It is reasonable to assume the quasiclassical equivalence 
in 3+1 case too. Therefore the main interest should be 
concentrated on the lowest states and the numerical pro- 
cedure we have used is mostly adequate namely for the 
lowest states. 

The resulting eigenvalues M 2 (K,y,a) were ordered in 
their magnitude and for each state separately the mini- 
mizing value of y, a was found. The orthogonality of the 
resulting wave functions for minimized parameters was 
then checked and appeared to be around few percents or 
better for all states except two. The latter belonged to 
almost degenerate states and the orthogonalization pro- 
cedure for them is made easily. 

The two-dimensional wave functions have been com- 
puted for minimized eigenvalues and formfactors and 
structure functions are computed according to the for- 
mulae of Section 4. 



VI. DISCUSSION OF RESULTS 

We have chosen six sets of quark mass parameters, 
including 4 sets of equal masses and 2 sets of unequal 
masses as shown in Table |. (all masses are in units 
of GeV). The values of parameters y,a,e are obtained 
by minimization procedure as explained in the previous 
section. The eigenstates obtained by our numerical pro- 
cedure, explained in Section 6, are listed in Table |n] and 
shown as Fig. [j]. 

Let us first discuss the light-cone spectrum. The 
Hamiltonian (|l|) depends explicitly on L z , whereas the 
Lorentz-invariance requires that masses do not depend 
on L z , but depend on L. Consequently one expects that 
the lowest state for L z = corresponds to L — 0, while 
the next state with L z — corresponds to L = 1 and 
must be degenerate in mass with the state L z = 1, L = 1. 
This is clearly seen in Table || and Fig.|l|. The next triplet 
of states with L z = 0,1,2 should correspond to L = 2 
Two of these states are found by us and minimized, while 
the third is not minimized since the procedure for this 
state was extremely time-consuming. As a result two 
states for L = 2 listed in Table || and Fig. [I] are degen- 
erated within 2% accuracy, while the third state is some 
12% up (but should come down after minimization). The 
same feature is roughly present for L = 3, 4 states where 



we have only one minimized state for each L. 

At this point one should take into account another de- 
generation - the dynamical one, and to this end discuss 
first center-of-mass spectrum, obtained from the cm. 
Hamiltonian j^] or by numerical diagonalization of the 
string equation of motion |J. We have computed the 
cm. spectrum using the routine from ]l5| . The cm. 
Hamiltonian for L = reduces to the so-called spin- 
less Salpeter equation [jl6| which was actually solved by 
us. The approximate form of the spectrum which can be 
easily computed also by the WKB method (jl7| can be 
represented as: 



(MW(L,N r ) 



■In a ( 2N r + -L + % ) (36) 

7T 2 



For L > one should take into account the string contri- 
bution, absent in spinless Salpeter equation but present 
in the Hamiltonian |j) and in 0. This correction found 
in § is: 



AM = - — 



16 2 L(L + 1) 



(M( ! 



o)) 3 



(37) 



The corrected masses, M(L,N r ) = M(°)+AM are shown 
in Table |l] and Fig.|l|. They are very close to the values 
computed in |^] . Now one can see that there is an ap- 
proximate degeneration in mass M = M^ ' + AM of 
states when one replaces one unit of N r by two units 
of L. This is seen in Table [n] and even better in 
fll8| |. The same type of approximate mass degeneration 
is seen in the light-cone spectrum - compare e.g. the 
states with (N r , L, L z ) = (1,0,0) ; M LC = 2.25 and 
(0, 2, 0) ; M LC = 2.28. 

In Fig.[j] this degeneration is visualized as the fact that 
masses appear on the vertical lines. This degeneration is 
a dynamical one and is a characteristic feature of non- 
relativistic oscillator. In our relativistic case it reveals a 
new string - like symmetry, typical for the QCD string 
spectrum |l6| ]. 

Let us now compare the light-cone and the center-of- 
mass spectrum. One could expect the coincidence up to 
an overall shift due to the different treatment of Z-graphs 
in two systems. In the cm. Hamiltonian these Z-graphs 
are presented but supposed to be unimportant on the 
grounds, that the backtracking of a quark in time nec- 
essarily brings about a folding in the string world sheet 
which costs a large amount of action and is therefore sup- 
pressed. The situation is different in the light cone - it 
is general belief that Z-graphs are absent here. One also 
expects that the Z-graphs (and the overall shift) should 
decrease if quark masses are increasing. 

The comparison of the spectra can be made from Ta- 
ble || and Fig.0. One can see indeed some overall shift 
down in the cm. spectrum by some 0.1 GeV and other- 
wise the masses coincide within the accuracy of compu- 
tation. This fact is highly nontrivial since two quantum 
Hamiltonians (the light-cone and center-of-mass ones) 



G 



are very different, they cannot be obtained from each 
other by a simple boost or other simple transformation. 
The light-cone Hamiltonian is rather complicated and 
it took the authors more than a year to get reasonable 
numerical results for it. 

We have also checked the quark mass dependence of 
the overall shift of spectra and proved that it drops 
sharply with the quark mass increasing see Fig. ^, sup- 
porting the idea, that the shift is due to different treat- 
ment of Z graphs (or self energy graphs). This fact also 
confirms, that the shift is not a consequence of some 
systematic errors of our procedure. In that case one 
should expect the decreasing of relative mass difference, 
i.e. SM/M with the quark mass increasing and not the 
decreasing of absolute difference SM, which has been ac- 
tually observed. 

We now turn to eigenfunctions. One expects in this 
case two types of excitation: the radial one leading to 
new nodes on p coordinate and similar to the 1+1 ex- 
cited states and the r± -excitation which causes nodes in 
the t-coordinate and associated with the orbital excita- 
tion. This is clearly seen in Fig. ^, where (a) refers to 
the ground state, (b) - to the orbital and (c) - to the ra- 
dial excitation. A more complicated example, combining 
both types of excitation is given in Fig. || (d) . 

As the next illustration we show in Fig. [| the case of 
two heavy quarks, demonstrating two types of excitation 
and also a new feature - the actual region of parameters 
is squeezed to a small region near p = 1/2. In Fig. | we 
demonstrate four states for the case of unequal masses - 
the "physical region" is shifted to one of the ends of the 
[0,1] interval. 

At this point it is important to find connection of our 
light-cone wave function to the nonrelativistic one, usu- 
ally defined in the cm.. In [|| it was demonstrated that 
this connection can be established only when both quarks 
are heavy, mi, m 2 ^ \fa. In this case y <§; 1 (as can be 
found by direct minimization of the Hamiltonian (^)) - 
the string transforms into the potential and loses its ma- 
terial and momentum contents. One can introduce as in 
[p[ the relative momentum p z and relative coordinate r z : 



Pz = (mi + m 2 ) x 



mi 



mi + m 2 



P+r_ 



mi + m 2 

and the M operator can be written as: 
1 



M « mi + m 2 + — (p 1 + p i) + C7 X r 
2m v 



(38) 



(39) 



(40) 



Hence one can write the momentum-space nonrela- 
tivistic wave function 2 ) directly through the light- 
cone variables: 



^(p 2 ) = * 



(mi + m 2 j 



m i 



mi + m 2 



This representation is valid in the large mass limit 
m, 3> \fo stated above and in addition near the center 
of the x-distribution, i.e. when \x — m\ / (m\ + m 2 )\ <C 1. 
The width of the peak in x variable is proportional to 
(mi + m 2 )~ 2 and is very narrow for heavy quarks. The 
form (|3^) is not correct for x at the ends of the interval, 
i.e. x = 0, 1 (remember that for large m^ the extremum 
value of y tends to zero and x = p changes in the inter- 
val [0,1]). Indeed the exact wave function as discussed 
in Section 4 vanishes linearly at x = 0, 1, while the r.h.s 
of (|4l]) stays nonzero. Moreover, the Jacobean J of the 
phase space d 3 p = Jdxd 2 p\_ is constant J — m\ + m 2 
and does not change this conclusion. 

The correspondence of the cm. and light-cone wave 
functions is lost if quark masses m, are of the order 
of -y/er or less. The physical reason is that the role 
of dynamics is now 100% important and the dynamics 
is different in different frames: the light-cone Hamilto- 
nian and wave functions are connected with those of the 
center-of-mass by a dynamical transformation, which in- 
cludes nonkinematical Poisson operators. Hence one can- 
not hope to obtain one wave function from another using 
only transformation including the kinematical (free) part 
of the operator M 2 . Moreover, a glance at the operator 
M 2 in ([!]) helps to realize that the separation of M 2 into 
the purely kinetic part (containing only momentum and 
masses) and purely potential one (containing only rel- 
ative coordinates) is impossible at all: e.g. the string 
enters through the term (P + r_) 2 and mixed terms like 
p±r± are present everywhere. This circumstance limits 
the use of simple recipes known in literature |?j], which 
connect the cm. and light-cone wave functions. In par- 
ticular, the ansatz for p z suggested in |jj and used for 
equal quark masses 



/ m 2 + fc 2 / _ 1 
x(l - x) \ ~ 2 



(42) 



coincides with our form (|39|) rigorously derived in |6|] only 
for large m, m^> \fa and x in the narrow region around 
x = 1/2 but differs at the tails of the wave function. For 
light quark masses, to < a and small p± the relation ( |42] ) 
yields incorrect results which can be seen in the almost 
constant behaviour in the interval [0,1] of the light-cone 
wave function, produced by the insertion of ([|2]) into the 
cm. wave function, typically 



to 2 + p\ 



1, 



(41) 



^(p z ) ~ exp(-aV 2 ) -> e M-a 2 [pl+ x{1 _^ (x--) 2 }) 

is insensitive to x for m,p± — > 0, whereas the exact light 
cone wave function is decreasing as x(l — x), see Fig. M, 
At the same time for heavy quark masses, m ^> %fa the 
two functions, one transformed by ( (42| ) from the cm. and 
another is a genuine light-cone wave-function solution of 
(|l3|), are very close to each other. This can be seen in 
Fig. @ for to = 1.4 GeV. 



7 



We now turn to the formfactor of computed states 1- 
4, see Table |, presented in Fig. ||. The main feature of 
the Fig. § is a very slow decrease of F(q 2 ) with q 2 which 
signals in particular a small radius of states. Indeed the 
values of V< r 2 > are too low (~ 0.338 fm for massless 
quarks). This fact is in agreement with the earlier cm. 
calculations of JTs) ] . 

Note that color Coulomb and spin interaction can only 
decrease radius. In the same Fig. || are also shown form- 
factor calculated in the cm. system for the equal mass 
values listed in Table | (cases 1-4). These formfactors 
can be calculated cither in terms of light-cone variables 
as in and |Tl"[ ], or simply using the nonrelativistic ex- 
pression 

F(q 2 )= J \^r)\ 2S ^j^d\ (43) 

since both integrals can be transformed one into another 
by a change of variables. One can see in Fig. || that 
the cm. formfactor is systematically below that of the 
light-cone. 

Finally we turn to the quark-distribution function 
q{p). It is computed through the light-cone wave- 
function tp(p,t) using ( p8| ) and shown in Fig. Q One 
can see the symmetric behaviour of q(p) with respect to 
reflection p — ► 1 — p. At the ends of the interval q(p) 
vanishes like p 2 and (1 — p) 2 , in the agreement with the 
1/q 3 behaviour of the formfactor at large q due to the 
Drell-Yan-West relations [||). 

Note the narrowing of the peak in q(p) in Fig. |9| for 
increasing quark masses. 

VII. CONCLUSION 

The present paper is the first in the planned series of 
papers devoted to the systematic study of nonperturba- 
tive contribution to formfactors, quark distributions and 
high-energy scattering amplitudes. 

The main physical idea of our approach is that the 
most part of nonperturbative dynamics in QCD is due 
to the QCD string, and the latter is described by the 
Nambu-Goto part of the Hamiltonian, which was writ- 
ten before in the cm. || as well as in the light-cone 
coordinates ||. 

Only valence part of the Fock's column was considered 
above in the paper, also for simplicity spins and perturba- 
tive gluon exchanges are neglected. To do the systematic 
comparison with experiment all these three simplifica- 
tions should be eliminated. Let us discuss their effect 
point by point. The higher Fock states are necessary 
to reproduce the Regge behaviour of q(x) ~ x ~ ap (®) a t 
small x (and at x — > 1 for quark distributions of hadrons 
in high-energy scattering). Here comes the first crucial 
point; to be answered in the second paper of this series, 
what is the QCD reggeon? 



In our method the higher Fock states, constituting 
the QCD reggeon, correspond to several gluons prop- 
agating in the nonperturbative background and there- 
fore confined to the excited Nambu-Goto surface [^o| . 
These states are in one-to-one correspondence with the 
excited Nambu-Goto string states. This is the picture at 
large distances; at small distances smaller than the vac- 
uum correlation length (width of the Nambu-Goto string) 
T g ~ 0.2 fm pi] , the string disappears and the usual 
perturbative gluon exchanges reappear. 

The effect of spin of light quarks is highly nontrivial 
||. It leads to the creation of the new vertex, which 
yields the constituent quark structure. Physically one 
may imagine this structure as being due to the light - 
quark walks around the end of the string. 

Having said all this, what is the lesson of the present 
work and of it possible development? 

The first lesson is that valence quark component can be 
successfully dynamically computed on the light cone; the 
string on the light cone is physically and mathematically 
well defined. The spectrum obtained on the light cone 
for the first time reasonably coincides with that of the 
cm. Hamiltonian for the string with quarks. Moreover, 
the nonperturbative wave function obtained directly on 
the light cone allows to calculate nonperturbative contri- 
butions to the formfactor and structure function. 

The second lesson is that the formfactor computed di- 
rectly on the light cone is close to of cm. for small q (see 
Fig. H), but is systematically above the cm. formfactor 
for larger q. This is not surprising, since in the light-cone 
formfactor there is a mechanism of the "redistribution" 
of the momentum q between the quarks, since q enters 
the light-cone formfactor ( p5| ) multiplied with (1 — x), so 
that the larger q, the smaller is (1 — a;) and the wave func- 
tion does not decrease too fast. Physically it means that 
at large q the configuration survives where the spectator 
quark gets as little momentum as possible so that it can 
be easily turned together with the active quark. This is 
exactly what is called the Feynman mechanism [jlj. 

The third lesson is that the minimal string plays only 
a passive role on the light cone, namely it participates in 
sharing of the total momentum and carries the part equal 
to < y >, but it docs not produce the x-distribution 
in structure function, which could simulate the gluonic 
structure function. The reason is that the string vari- 
able - the einbein field v - is quasiclassical and has no 
dispersion. 

The value of < y > computed according to (pl|) de- 
pends on quark masses and is listed in Table. Hit is 
resonable that < y > = 0.22 is smaller than the experi- 
mental value of overall gluon momentum, 0.55, since in 
our picture the difference should be filled in by higher 
Fock components. 

The fourth lesson comes from the comparison of the 
computed quark distribution, Fig. ^, with the experi- 
mental data for the pionic structure function p2| . Be- 
haviour of q(p) at small p and small (1 — p) is symmetric 
in Fig. ^, while in reality q(x) should rise at small x like 
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x - a p(u) ^ x~ ' 5 (we neglect at this point the difference 
between x and p, which is due to < y >). This peak 
at the small x should be filled in by the contribution of 
higher Fock components, containing additional gluons on 
the string, as was discussed above. The behaviour of q(p) 
at p — 1, which is calculated to be (1 — p) 2 will be also 
changed into (1 — p) due to gluon exchanges, which ac- 
count for the formfactor asymptotic 1/q 2 at large q, and 
the Drell-Yan-West duality ensure the (1 — p) behaviour 
around p = 1. 

Finally, the formfactor calculated above in the paper, 
Fig. ^ shows too little radius of the "pion" < r 2 
(0.34 fro) 2 as compared with the experimental one < 
r 2 >~ (0.67 fm) 2 . This fact is in qualitative agreement 
with other calculations, where the cm. wave-function 
was used [^9|, and some authors assumed as in that 
quarks have "internal" structure and their own radius 
which should be added to the "body radius" to repro- 
duce the experimental value. This fact of small body 
radius seems to be a necessary consequence of the simple 
string + point-like quarks picture, and probably cannot 
be cured by the higher Fock components. 

The calculation of the quark structure as produced by 
the spin and chiral effects is thus an interesting and fun- 
damental problem which will be discussed in another pa- 
per of this series. 



cm. Hamiltonian the role of T is played by the center- 
of-mass euclidean time coordinate T = (x^ + x^)/2 and 
the hypersurface is a hyperplane £4 = X4 — const. 
With the notations for the vectors (2^,6^ 

ab = a^b^ — aibi — ao&o — <x_l&_l + a + 6_ + a_& + , 
a 3 ± a 

one can define the hypersurface ^2 through the qq coor- 
dinates z^.z^ as 

and the kinetic part of the action A 

A = K + K + aS m in, 

has the form 

K + K = - A J^ zl(r)dr+ 

1 f s 2 f 3 f* 

+— / !„(f)cif+ / m 2 dr + / rn^dr = 

4 Jo Jo Jo 
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dz 4 



H ( z 1 + 2z-)-\ -H - 



where we have defined 

2 Ml (z + ) = ^± ; 2n 2 (z + ) = ^± 

The minimal surface S m in is formed by connecting 
Zfj,(z+) and z^{z+) with the same value of the evolution 
parameter z + , i.e. 



APPENDIX A: 

This appendix is based on the material of Ref. [|| 
Given a qq Green's function in the coordinate space 
G(xx;yy), where xx(yy) are final (initial) 4-coordinates 
of quark and antiquark, one can define the Hamiltonian 
H through the equation (in the euclidean space-time) 



dG 
df 



= -HG 



(A.l) 



where T is an evolution parameter corresponding to some 
choice of a hypersurface S. In a particular case of the 



rr \ dz + ( d{3[w 2 W 12 - (W) 2 ] 1/2 











where 



w^(z+;f3) = z^(z + )(3 + z^(z + )(l - (3) 



and dot and prime denote derivatives in z+ and (3 respec- 
tively throughout this Appendix. 

We now introduce " center-of-mass" and relative coor- 
dinates, 

Rp xz^ -\- (1 x)z^ , T ^ — z^ z ^ 







where the variable x is defined from the requirement that 
the term r±R± should be absent in the action. This 
yields: 

Hi + J vfidfi 



Mi + H2 + f vdj3 
Then for the Green function one obtains: 

G(xx; yy) = J Dfii(z + ) D[i 2 (z + )Di> DR^ Dr^ e" 



where the action A 

if mi 
- + — 

Mi M2 



A = - I dz + { -i + + ai {R\ + 2R-) + a 3 rj 



and r_. As it was shown, canonically conjugated mo- 
mentum to the i?_ is P+ = oi. Situation with r_ is 
more subtle. In order to clarify the situation let us start 
with the general form of the qq Green's function in the 
Feynman-Schwinger formalism 

G(x,y) = J ds dsDzDz l - K - R < W{C) > 

to impose boundary conditions, one can rewrite DzDz 
using discretization 

= z{n) - z(n - 1), = z(n) - z(n - 1) 
DzDz — Y[ d(, n d£ n > dp dp'x 



, o 2 (r_+R ± r ± + (<[3>-x)r ± r ± ) 2 
+ / — dp rj_ 



(A. 



(r±r ± ) 2 / v{fi- < (3 >) 2 d(3 



The following notation was used: 



ai= m + j-i2 + / v([3)df3 
Jo 



a 3 = m (1 - x) 2 + il 2 x 2 + / v{(3)((}- x) 2 dP 

Jo 

Integration over DR^ leads to an important constraint: 
a x =P+ 

Furthermore we go over into the minkowskian space, 
which means that 



■ M ■ M 

Hi — > -i^i , v — > -iu 



M 



ai^-iaf, A^~iA 

For the minkowskian action we obtain (omitting from 
now on the superscript M everywhere) 



1 



A M = ~ dz + -^i-^ + asri 



Mi 



a 2 d(3 2 



v ai 



rj_ (/ii+M2)H 



[r_ + (< (3 > ~x)r ± r ± ] 



and using (4) one easy obtains the Hamiltonian (Q) which 
opened the main text. 

To complete the Hamiltonian formulation of our prob- 
lem we define canonical momenta for the coordinates R- 



x exp [ip(Y^ i n + X-Y)+ ip'{Y, In' + X - Y) } 
One can introduce the total and relative momenta 



P=p + p', q = 



p-p 



and R — AR n /e; Ne = T, one has 

£n*^n £n(I %n) AR n 
£n £n — AVn 

Expressing through AR n ,Ar n and going over to 

the momentum representation of G one obtains 

G(P) = J dq]J dAR n dAr n , x 
xexpl iA + iP^2 AR n + i^2(^P(l - 2x) +q)Ar n \ 

\ n n ) 

where A = J Q T drC = J Q T dr(K + K — aS m i n ) One can 
now introduce the Hamiltonian form of the path integral 
via 



Dxe 



fCdr = 



DxDpe 



and rewrite the first two exponents as 

exp^i J PiR.dr + i J [^Pi(l - 2x(t)) + qi ]hdT 

The term proportional qi disappears because of boundary 
conditions r M (0) = r^(T) — 0, and one obtains 

P+ = \p+{1-2x), [p+,r_] = -t 

One can rewrite this in the form 
d 

P+r_ = i—, [P+r_, x] = i 
ox 
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APPENDIX B: 



The aim of this appendix is to derive the Eq. (|l^) from 
Eq.(nl). The first thing to do is to express the Hamilto- 



nian fllf) as 



H 



M 2 
2P+ 



(B.l) 



Then using definition (|3j) we obtain the following form 
for M 2 : 



M 2 = 



mi TO, Lz 
1 1 — 

x-y </3> 1 -x -y (l- < (3 >) ~cr\ 



+ 



1 



x-y a <f3> l-x-y {l- < P>) 



+ / ^r 2 ± df3 



1 - Vo r 2 



(B.2) 



here 



Quantization of the above expression is done according 
to the canonical commutation relations: 

{p k ± ,r j ,}=-i6 k i 



{x,(P+r_)} = -z 

We are looking for the wave function of the problem given 
in the mixed coordinate - momentum representation ip = 
ip(p, t), so one has to substitute into (B.3) the operators: 



(P+r-)=i 



1 



d 



i - m J dp 



p'i 



d 



The important point is the operators ordering. We use 
the Weil ordering rule, i.e. 

AB — ► ^(AB + BA) 

for any noncommuting operators A and B. Let us notice 
that y explicitly depends on t according to (||) and hence 
should also be differentiated during the ordering proce- 
dure. The final result for the operator M 2 may be found 
by the straightforward calculations, it is 

rtt a 92 a 92 a 92 

M = Ai —J + A 2 —J + A 3 —— + 
dr dp 2 - otdp 



and 



c=(x-y o <0 >)(! -x) 2 + 



[1 - x - yo(l- < >)]x 2 + y{(3- xfdfi 

Jo 



7 = ~r^~~ r ~ (<P> 



i -yo 

As it was mentioned above, it is more convenient to ex- 
press M 2 through a new variable p: 



1 



P = 



— (a; -2/0 <P>) 



1 - yo 

Substituting this variable into (B.2) one obtains: 
M 2 = 1 fm 2 : m 2 L 2 



1 - yo L p i-p ct 



l l 



p i - p 



[p ± r ± +7(P+r_)] 2 



y t 



where the notation t = r\ was used. 



(B.3) 



A 9 A 9 

+A 4 — + A 5 — 
at op 



where the coefficients A, 



A 1 (p,t) = - 



At 



A, 



l + at 



(B.4) 



p(l-p) [l + (a-y)t] 



MM 



y(l + at) 

[i + (a - y)t)] 3 



yt ^.^f + l + at 
Pi 1 - P) 



M{p,t) = -4yt 



(1 + at) (p- </3>) 
[l + (a-y)t] 2 p(l-p) 



x <^ 4 



M(p,t) 



yt 



l 



p(i-p) 



l + at 



[1 + (a - y)t] 2 1 + (a - y)t 



2yt 



{l + at) (p 2 -2p</3> + <{3>) 



[1 + {a - y)t] 



P(l-P) 



A 5 (p,t) = -2y 



[1 + (a + y)t] (p- </3>) 
[l+(a-„)t]3 p(l-p) 
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2 (1 + at) 
V [l + (a-y)<F X 



(p-<f3>)[(p(l-2<f3>)+<0>] 



W-pW 



Mp,t) = 



1 + at ( m\ mj \ 1 + at 



l+(a-y)t \ p 1- p) y 



L 2 z (l + at) 2 y 2 t(l + at) 



t a 6 [1 + (a - y)t] 3 



[p 3 (l - 2 < >) - 3 < > 2 p(l - p)+ < (3 > 2 



[P(l-P)] 3 



y [3 -(a- y)t] 
'p(l-p) [l + {a-y)tf 



{p 2 - 2 < (3 > p+ < (3 >) [1 + {a + y)t] 



W-pW 



[l + (a- y)tf 



where: 



a 6 = yt[l + {a - y)t]{p- < (3 >) 2 + 



+ [1 + {a - y)t](l + at)p(l - p) + yt{l + at)j 
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FIG. 1. The Chew-Frautschi plot with masses computed 
via the light-cone (circles, squares and triangle for L z =0,1,2 
respectively) and the cm. Hamiltonian (stars). The system- 
atic overall mass shift is seen as a divergence of straight lines 
passing through circles and stars. The states with high L or 
high AT r (daughter trajectories) are numerically less reliable 
and not shown. 

FIG. 2. The systematic mass shift between mass eigen- 
values of the light-cone and cm. Hamiltonians versus quark 
mass. 

FIG. 3. The 3d plots of wave functions of the four lowest 
states of light-cone Hamiltonian for zero quark masses. Co- 
ordinates on horizontal plane are < p < 1 , < t < 15 (in 
units of o --1 ). 

FIG. 4. The same as in Fig. 3 but for heavy quark masses, 
mi — mi = 5 GeV. 
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FIG. 5. The same as in Fig. 3 but for unequal quark 
masses, mi = 5, ni2 = GeV. 

FIG. 6. The 3d plots of the ground-state wave func- 
tions computed via the light -cone Hamiltonian (up- 
per part) and via the cm. Hamiltonian, with the standard 
substitution (E^) (lower part) for zero quark masses. 



FIG. 7. The same as in Fig. 6 but for heavy quark masses, 
mi = mi = 1.4 GeV. 

FIG. 8. Formfactors calculated with light-cone 

wave-functions for the cases 1-4 of Table | (solid lines) and 
with the cm. wave-functions (lines with stars). 

FIG. 9. The quark-distribution function q(p) computed 
with light-cone wave-functions for the cases 1-4 of Table fl 



TABLE I. Quark masses and minimizing values of vari- 
ational parameters for six cases, computed in the paper to- 
gether with < r 2 > and < y > for each case. 





m qi 


m 92 


y 


a 


e 




(y) 




Case 1 


0.00 


0.00 


0.40 


0.75 


1.0 


0.50 


0.223 


0.338 


Case 2 


0.25 


0.25 


0.35 


0.75 


1.0 


0.50 


0.190 


0.329 


Case 3 


1.40 


1.40 


0.20 


0.75 


1.7 


0.50 


0.081 


0.249 


Case 4 


5.00 


5.00 


0.10 


1.50 


3.5 


0.50 


0.017 


0.1675 


Case 5 


1.40 


0.00 


0.20 


0.50 


1.1 


0.10 






Case 6 


5.00 


0.00 


0.07 


0.25 


1.5 


0.13 







TABLE II. The light-cone mass eigenvalues (mass in L-C 
- the fourth column), and the cm. mass eigenvalues without 
(mass in C-M - the fifth column) and with string correction 
(M (0) - AM in CM - the sixth column) The first 3 columns 
contain the quantum numbers assignment for the given state, 
and the last two columns - the minimizing values of varia- 
tional parameters y and a. All eigenvalues are computed for 
quark masses equal to 0.12 GeV. The masses marked by as- 
terix are minimized with parametres y, a listed in the last 
two columns. 





L 


L z 


Mass 






M (0) - AM 


y 


Q 








m Li— \u 




m o— ivi 


m o— ivi 






U 


A 
U 


n 
U 


1.0 1 U4 




1.4uU4 


1.4DU4 


n /in 

U.4U 


n on 





1 





2.0063 


* 


1.9236 


1.8637 


0.15 


0.20 





1 


1 


1.9926 


* 






0.40 


0.80 





2 





2.2835 


* 


2.2990 


2.1937 


0.20 


0.30 





2 


1 


2.6138 








0.40 


0.80 


U 


n 

z 


o 
z 


Z.oZZZ 


* 






n /in 
U.4U 


U.oU 


U 




U 


O RK*7 A 

2. DO 1 4 


* 






U.1U 


n on 
U.oU 





3 


1 


3.1665 








0.40 


0.80 





3 


2 


2.8157 








0.40 


0.80 





4 





3.0675 


* 






0.15 


0.20 





4 


1 


3.3175 








0.40 


0.80 





4 


2 


3.3127 








0.40 


0.80 


1 








2.2567 




2.1483 


2.1483 


0.20 


0.40 


1 


1 





2.6477 


* 


2.4728 


2.4446 


0.10 


0.15 


1 


1 


1 


2.5909 








0.40 


0.80 


1 


2 


1 


3.0787 








0.40 


0.80 


1 


2 


2 


2.8879 








0.40 


0.80 


2 








2.9298 




2.6707 


2.6707 


0.10 


0.20 
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